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SUMMARY 


The  effort  was  initiated  by  a mathematical  treatment  of  the  laws  of  mass 
flow  and  energy  conservation  and  the  momentum  theorem  as  they  apply  to  a model 
representing  an  ejector  in  which  the  processes  of  ingestion  and  mixing  occur 
sequentially  in  separate  and  distinct  regions  of  the  duct.  Injection  is  assumed 
to  occur  at  the  start  of  the  mixing  section,  downstream  of  the  ejector  inlet. 

Mixing  occurs  in  a constant  cross-section  duct  after  which  the  flow  is  exhausted 
in  an  outlet  diffuser  or  nozzle  to  ambient  pressure.  Jet  diffusion  is  applied 
if  practical. 

Injected  and  induced  flows  are  assumed  to  consist  of  compressible,  ideal 
fluids  having  arbitrary  properties  consistent  with  the  perfect  gas  law.  The 
losses  due  to  inlet  blockage,  skin  friction  and  incomplete  thermal  and  kinetic 
mixing  were  teUten  into  consideration  in  the  formulation  of  the  mathematical  system 
of  equations.  In  addition  provision  was  made  to  assess  the  influence  of  boundary 
layer  ingestion  into  the  ejector. 

Solutions  of  this  set  of  equations  were  based  upon  assumed  values  of  the  area 

ratio  a (=  X-/a_) , the  environmental  characteristics  (U  , p„,  T„  and  M ) , the 
OO  2 ^ 0 *00'  00  00  ' 

properties  of  the  injected  gas  (m  , T and  P ) , the  magnitude  of  the  loss  coef- 

_ P op 

ficients  (C, .,  n.»  n„,  C,,  As),  the  non-uniformity  due  to  incomplete  mixing  and 
Gi  1 N r 

the  pressure  of  the  induced  flow  at  Station  1 (Pj^) . With  these  given  quantities, 

the  chciracteristics  of  the  flow  at  any  station  of  the  ejector  and  the  outlet  area 

ratio  5 (=  X^/X^)  were  determined  from  the  set  of  equations.  (The  pressure,  p^^, 

can  be  calulated  by  the  same  set  of  equations,  if  6 is  specified) . The  performance 

parameters  (({i,  n , R ,11  » H..)  could  then  be  evaluated  in  terms  of  the  flow  charac- 
p p o t 

teristics  for  emy  given  set  of  input  quantities. 

When  the  density  of  the  fluid  does  not  vary  with  pressure  or  temperature 
(as  in  the  case  of  an  incompressible  fluid)  the  ejector  flow  equations  have  a 
single  valid  solution.  The  functional  dependence  of  density  upon  pressure  and 
temperature  in  a compressible  fluid  however,  results  in  a situation  in  which 
two  solutions  can  represent  physically  realistic  flow  properties. 

The  first  solution  (under  which  the  flow,  after  complete  mixing,  is  always 
subsonic)  provides  information  reg^u:ding  the  flow  pattern,  performance  and  ejector 
geometry  corresponding  to  any  given  flight  or  operational  condition  and  injected 
gas  characteristics.  Under  this  solution  the  resulting  thermodynamic  properties 


1 


of  the  mixed  flow  always  are  consistent  with  the  mass  flow  and  energy  conser- 
vation laws  and  with  the  momentum  theorem  and  the  Second  Law  of  Thermodynamics, 

An  optimization  procedure  in  which  the  value  of  the  Mach  Number  of  the  induced 
flow  at  the  start  of  mixing  (Mj^)  was  systematically  varied  over  a wide  range, 
indicated  that  the  ejector  performance  was  a maximum  at  one  value  of  which 
depends  upon  the  values  of  the  other  input  parameters.  Results  of  this  numerical 
analysis  indicated  that  the  optimal  value  of  under  this  first  solution  is 
always  less  th^m  1.0,  but  varied  from  values  close  to  1.0  to  low  subsonic  values, 
depending  upon  the  flight  Mach  Number  emd  the  characteristics  of  the  injected  gas. 
Ejector  performance  at  this  design  point  decreased  from  high  values  of  thrust 
augmentation  under  stationary  conditions  to  values  close  to  1.0  at  the  high 
subsonic  flight  Mach  Numbers,  when  the  injected,  energized  gas  was  unheated. 

With  the  injection  of  hot  gas  without  mechanical  compression,  the  ejector 
performance  obtained  from  ejectors  designed  under  this  first  solution  with 
subsonic  mixing  remained  very  high  over  the  entire  range  of  practical  flight 
Mach  No's.  Supersonic  mixing  provided  high  performance  only  at  supersonic 
flight  speeds.  The  injection  of  gases  having  high  thermal  and  mechanical 
energy,  as  does  the  efflux  from  a turbo-jet,  resulted  in  a decrease  of 
ejector  performance  similar  to  that  with  unheated  injected  gas  as  the  flight 
Mach  Number  varied  from  zero  to  about  0. 5,  Further  increasing  the  flight  Mach 
Number  produced  improved  performance  throughout  the  remainder  of  the  practical 
flight  speed  rcmge  with  ejectors  designed  under  this  first  solution,  with  high 
temperature  injected  gas.  Thus  it  is  evident  that  at  high  flight  speeds, 
ejectors  operate  effectively  with  hot  gas  injection  under  the  first  solution. 

The  second  solution  (under  which  the  flow  after  complete  mixing  is  always 
supersonic)  provides  information  regarding  the  flow  pattern,  ejector  performance 
and  geometry  corresponding  to  any  given  flight  or  operational  condition  and 
injected  gas  characteristics.  Variation  of  over  a wide  range  indicates  that 
the  total  entropy  of  the  completely  mixed  flow  may  be  less  than,  equal  to  or 
greater  than  the  sum  of  the  total  entropies  of  the  primary  and  induced  flows 
prior  to  mixing.  Low  subsonic  values  of  M^^  generally  produce  negative  changes 
of  total  entropy  which  increase  as  is  increased,  becoming  greater  than  zero 
prior  to  M^^  = 1.0.  Thus  at  some  subsonic  value  of  the  entropy  change  (S2“ 
is  zero.  This  point  is  chosen  as  one  of  two  characteristic  points  under  the 
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second  solution.  Continuing  to  increase  results  in  positive  entropy  ch2mges 
and  in  general,  a functional  dependence  of  thrust  agumentation  ((())  upon  which 
has  a local  maximum  at  a value  of  which  may  be  slightly  less  th£m  1.0,  but 
which,  in  general,  is  greater  than  1.0.  This  maximum  ((i  is  considered  to  be  the 
second  characteristic  point  under  the  second  solution. 

In  general  there  are  three  distinct  characteristic  points,  one  which  results 

from  the  use  of  the  first  solution  and  two  of  which  result  from  the  use  of  the 

second  solution,  all  of  which  are  consistent  with  the  Second  Law  of  Thermodynamics. 

The  performance  and  the  geometry  of  the  ejector  are  determined  by  the  flow  charac- 
teristics at  the  design  point  or  in  other  words  <p,  Hp  ^ j » ^ j > Hq  gj 

determined  by  the  value  of  M, (or  p, ) , the  input  quantities  (a  , M , AP/p  , AT/T  ) 

X.  00  00  OO  00 

and  the  calculated  quantities  (a,  6 , X , X,,  V /V  ). 

J poo  pi 

Examples  illustrating  the  dependence  of  (}>  upon  are  presented  for  both 
solutions  at  two  operational  conditions  (M^  = 0.3  and  = 2.0)  using  injected 
gas  characteristics  corresponding  to  typical  effluent  gases  from  a fanjet  at 
= 0.3  and  from  a ramjet  at  = 2.0,  to  illustrate  the  shape  of  the  curves, 
the  location  of  the  three  characteristic  points  eind  the  criticality  of  the  proper 
choice  of  or  ejector  configuration. 

Ejector  performance  over  the  entire  practical  range  of  flight  speeds  is  then  'c 

presented  for  fanjet,  ramjet,  turbojet  and  rocket  type  gas  generators  at  each  of  ' 

the  three  characteristic  points,  to  illustrate  the  influence  of  combination 
thrusters  on  the  performance  of  each  type  of  gas  generator  acting  as  a free  jet. 

The  use  of  a jet-diffuser  modifies  the  set  of  equations  since  the  flow  does 
not  return  to  ambient  pressure  at  Station  3 as  in  the  case  when  diffusion  is  ‘ 

accomplished  by  solid  surfaces.  The  diffuser  jet  itself  requires  the  injection 
of  energized  fluid  whose  momentum  must  be  considered  in  the  evaluation  of  ejector 
performance.  The  relationship  among  the  diffuser  jet  slot  thickness,  angle  of 
departure  from  the  solid  surface,  the  effective  area  ratio  and  the  momentum 
requirement  are  presented  for  use  in  evaluating  the  performance  of  this  type  of 
ejector. 

The  real  fluid  effects  and  the  influence  of  incomplete  mixing  are  then 
discussed  and  compared  with  existing  experiments  performed  under  stationary 
conditions.  The  comparisons,  using  realistic  values  for  the  loss  coefficients 
indicate  agreement,  proving  the  validity  of  the  method  and  the  relative 
unimportance  of  the  requirement  for  complete  mixing. 
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SECTION  I 


INTRODUCTION 


Basically  ejectors  are  suitably  shaped  ducts  which  are  designed  to  capture 
ambient  fluid  and  cause  it  to  mix  with  a stream  of  energized  fluid.  This  concept 
was  first  reduced  to  practice  in  the  form  of  a pump  which  caused  a large  mass 
flow  of  essentially  stationary  fluid  to  flow  against  a small  pressure  head  while 
being  energized  by  a small  mass  flow  of  rapidly  moving  fluid.  Recognition  of 
the  fact  that  the  ejector  process  could  provide  an  increase  of  momentum  flux 
compared  to  that  of  the  energized  stream,  stimulated  the  development  of  ejector 
thrust  augmenters.  Early  theoretical  analyses  based  upon  the  theory  of  flow  of 
incompressible  fluids  appeared  to  indicate  a large  potential  for  augmenting  the 
thrust  of  a fluid  stream  when  the  ambient  fluid  was  induced  to  flow  through  the 
ejector  from  a stationary  undisturbed  state  or  when  the  ejector  was  at  rest  with 
respect  to  the  undisturbed  fluid.  Motion  of  the  ejector  in  its  thrust  direction 
appeared  to  result  in  a rapid  decrease  of  net  thrust  with  increasing  ejector 
velocity.  This  observation,  based  upon  some  limited  investigations  which  did 
not  attempt  to  optimize  the  thrust  augmentation  by  proper  selection  of  ejector 
geometry  and  injected  fluid  characteristics,  discouraged  further  attempts  to 
utilize  ejectors  as  thrust  augmenters  during  translational  motion.  In  addition, 
since  the  entrainment  of  ambient  fluid  was  known  to  depend  upon  mixing  of  the 
energized  stream  with  the  ambient  fluid,  it  was  assumed  that  complete  mixing 
was  essential  regardless  of  the  cost.  This  assumption  led  to  the  design  of 
ejectors  having  large  cross-sections  to  capture  large  amounts  of  ambient  fluid, 
and  multiple  injection  jets  of  energized  fluid  to  accelerate  mixing.  Further 
it  was  recognized  that  by  incorporating  a diffuser  at  the  end  of  the  mixing 
process,  it  was  possible  to  reduce  the  throat  pressure  and  to  convert  lost 
kinetic  energy  into  useful  energy  thereby  increasing  the  flow  of  ambient  fluid 
per  unit  area  and  improving  the  propulsive  efficiency.  The  use  of  diffusers 
however  resulted  in  further  increases  in  ejector  length  and  diameter  as  well  as 
the  momentum  loss  due  to  skin  friction  and  separation.  Furthermore,  if  flow 
separation  occurs  in  the  diffuser,  the  diffuser  area  ratio  is  effectively 
reduced.  As  a matter  of  fact,  without  boundary  layer  control,  effective 


diffuser  area  ratios  in  excess  of  2 do  not  exist  in  ejector  technology.  Thus 
ejector  thrusters  were  considered  to  be  impractical  on  the  basis  of  three  concepts 
derived  primarily  from  incompressible  flow  theory. 

1)  Decrease  of  augmenting  capability  with  increasing  translational 
velocity. 

2)  Long  length  required  for  complete  mixing. 

3)  Large  transverse  and  longitudinal  dimensions  required  for  effective 
diffusion. 

Recent  advances  in  diffuser  design  have  resulted  in  the  development  and 
testing  of  a jet-diffuser  (Reference  1)  vdiich  has  the  ability  to  achieve  large 
diffuser  area  ratios  with  small  lengths  and  lateral  dimensions  compared  to  those 
of  solid  diffusers.  This  concept  involves  the  use  of  a high  velocity  thin  jet  of 
fluid,  injected  at  the  wall  of  the  solid  diffuser  surface  to  avoid  separation  and 
to  provide  a curved  jet  sheet  downstream  of  the  solid  surface  to  sustain  a pressure 
gradient  and  thus  extend  the  region  in  which  the  flow  can  diffuse  and  mix. 
Experiments  have  verified  that  solid  surfaces  downstream  of  the  diffuser  jet  can 
diverge  to  half  angles  (B)  at  least  as  large  as  60  degrees  without  separation. 

Investigations  of  the  optimized  performance  of  solid  diffuser  ejectors  in 
high  speed  motion  in  their  thrust  direction  in  an  incompressible  fluid  have  been 
carried  out  recently  by  FDRC  (Reference  2)  where  it  has  been  shown  that  the  thrust 
augmentation  does  not  necessarily  decrease  with  increasing  speed.  In  fact  it  has 
been  shown  that  in  em  incompressible  fluid,  thrust  augmentation  increases  with 
decreasing  values  of  q/AP,  the  ratio  of  the  dynamic  pressure  to  the  pump  head. 

Thus  regardless  of  speed,  the  thrust  augmentation  remains  constcint  when  q/AP  is 
consteuit.  In  other  words  if  pump  heads  (AP)  increases  in  proportion  to  the 
dynamic  pressure  (q)  , the  thrust  augmentation  (<(>)  remains  constant,  provided 
other  factors  such  as  the  inlet  emd  diffuser  area  ratio  are  fixed  and  cavitation 
does  not  occur.  This  work  was  extended  to  include  the  jet-diffuser  ejector  and 
is  reported  in  Reference  3,  where  it  is  shown  that  the  thrust  augmentation  depends 
upon  the  ratio  q/Ap  as  in  the  case  of  a solid  diffuser  ejector. 

With  the  background  of  the  incompressible  flow  analyses,  the  effort  reported 
in  this  document  was  undertaken  to  determine  if  the  performance  of  ejectors  in  a 
con^ressible  medium  could  be  optimized  in  terms  of  the  injected  fluid  character- 
istics so  that  high  performance  could  be  achieved  at  high  speeds  with  selected 
types  of  gas  generators. 
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SECTION  II 
TECHNICAL  DISCUSSION 


1.  Analytical  Model 

The  analysis  of  ejector  flow  properties  and  the  evaluation  of  performance 
in  a compressible  fluid  environment  has  been  carried  out  under  the  assianptions 
that 

1)  All  fluids  involved  in  the  processes  are  perfect,  compressible  gases. 

2)  Injected  fluids  are  at  arbitrary  pressure  and  temperature. 

3)  Induced  flow  is  at  an  arbitrary  velocity  and  stagnation  pressure. 

4)  Influences  of  inlet  blockage  or  shock  losses,  skin  friction  and 
diffuser  seperation  are  included  in  the  analysis. 

5)  Mixing  occurs  in  a region  of  constant  cross-section. 

6)  Jet  Diffusion  is  used  where  practical. 

The  schematic  arrangement  of  this  generalized  ejector  configuration  with 
alternative  inlets  and  outlets  and  the  pertinent  station  designations  are 
illustrated  on  Figure  1. 


pi)  (Primary  Flow) 


(Induced  Flow) 


In  general  the  equations  of  state,  the  laws  of  energy  cind  mass  flow 
conservation,  the  momentum  theorem,  and  the  Second  Law  of  Thermodynamics 
have  been  imposed  to  obtain  valid  performance  estimates  and  flow  field 
configurations. 


2.  Analysis 

The  thrust  efficiency  of  the  primary  nozzle  (n^j)  is  represented  by  a 
reduction  of  the  ideal,  iaentropic  velocity.  Thus 

V = n V ' 
p“  N p“ 

where  the  primed  velocity  is  the  velocity  which  would  result  from  an 
isen tropic  expansion  from  the  ideal  stagnation  pressure  to  the  ambient 
pressure  (p^^)  . Thus  since 


(1) 


(2) 


and 


vVvop[^-(y”] 


(3) 


op' 


where  n = (y-D/Y 

P is  the  ideal  lossless  reservior  pressure  auid 
or 

P is  the  stagnation  pressure  of  the  primary  nozzle  (with  loss) . 
op 


The  thrust  efficiency  can  then  be  expressed  as 


The  Inlet  to  the  ejector  is  assumed  to  have  a velocity  CJ^  and  a stag- 
nation pressure  ^Ich  differ  from  the  freestream  velocity  and  the 
isentroplc  stagnation  pressure  P^,  eui  a result  of  the  boundary  layer  or  wave 
formation  associated  with  the  vehicle  drag.  However  the  stagnation  temper- 
ature is  Identical  to  that  of  the  freestream  (T^^  = T^) . Assuming  the  velocity 

ratio  at  the  ejector  inlet  Is  U Al„»  it  can  be  shown  that 

0 00 


(6) 


(7) 


The  skin  friction  at  the  surfaces  of  the  inlet  of  the  ejector  results  in 
a further  loss  which  is  represented  by  an  inlet  loss  coefficient  which 
is  defined  as 

-oi  ■ <V2'“i  =ai  ■ 


which  reduces  to  the  convenient  form 
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This  concept  is  useful  when  viscous  drag  is  dominant.  However  when  wave 
drag  is  dominant  the  inlet  loss  is  more  conveniently  expressed  as  an  inlet 
recovery  factor  (ri^^)  where 


P 

oe 


(9) 
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The  mixing  process  is  analyzed  under  the  assumptions  that  the  process 
occurs  in  a constant  area  duct  between  Stations  1 emd  2,  from  known,  uniformly 
distributed,  arbitrary  properties  at  Station  1 to  a final  state  at  Station  2 
where  the  mixing  is  only  partially  completed.  The  mixing  process  is  considered 
as  being  comprised  of  two  processes}  the  thermal  mixing  in  which  the  density 
distribution  at  Station  2 is  described  as 


and  the  kinetic  mixing  in  which  the  velocity  distribution  at  Station  2 is 
described  as 


where  the  primed  terms  are  assumed  to  be  small  compared  to  the  average  guemti 
ties.  Pressure  at  Station  2 (p^)  is  assumed  to  be  uniformly  distributed  in  a 
transverse  direction,  emd 


Under  these  conditions,  the  law  of  energy  conservation  is  expressed  as 


in  view  of  the  equation  of  state 


r t 


The  momentum  theorem  including  the  influence  of  a frictional  force  {f  ) due 


to  skin  friction  between  Stations  1 and  2 is  expressed  as 

P-y  * ! Pt  “ Pi  " V + m.u 

22  '22-'  m *^12  ppl  il 

^2 

where  f^  = Force  due  to  skin  friction  in  the  mixing  duct  (2-dimensional) . 

If  Z = Length  of  mixing  duct 
rn 

and  Coefficient  of  skin  friction  in  the  mixing  duct. 


(17) 


then 


and  Equation  17  is  more  conveniently  written  in  a dimensionless  form  as 
(P2  - Pi)« 


~ ^ ^ ^1  - srh:  ^ ‘^2  ^ pp  <'^2  - S 


- ,,2  - - m V , ^ '^2  ''2'  '”2  ' “2'  ''fm 

'pl'^pl  P Ppl'l^l 


P2  "2  ^2 


(18) 


To  evaluate  the  pressure  terms  in  the  momentum  equation,  consider  the  total 
enthalpy  at  Station  2. 

2 

(23) 


Y ^2  . ^2 


Y - 1 P,  2 Vo2 


which  when  mass  averaged  over  the  cross-section,  under  the  assumed  conditions 
becomes 


1 Pj  / ^ • i / 1^2  * . c / T^2« 


(24) 


X2  2 

or  since  dm  = P^U^dy,  this  becomes 


7^  P2  / U^dy  J / (P2  + p')  (U^  t up  dy  = Cp  / T^^dm  (25) 

^2  \ ^2 
and  upon  integration,  neglecting  third  or  higher  order  terms,  this  yields. 
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(26) 


Using  Equation  20  or  21  to  eliminate  the  mass  flow  terms  results  in  the 
expression 
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which  can  also  be  expressed  as 


or 
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where  J 
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The  remaining  pressure  term  in  Equation  22  can  be  expressed  conveniently 
using  the  equation  of  state  and  the  enthalpy  relationship  as 


Pl“ 
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Therefore 
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Using  the  above  expression.  Equation  22  reduces  to  the  form, 

,2 


A X 


B X^  + C = 0 


(33) 


This  is  a quadratic  equation  whose  solution  is 
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noting  that, 
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In  the  absence  of  slcin  friction  and  with  complete  mixing,  the  two  solutions 
of  Equation  34  provide  Mach  numbers  (M^^  and  M^j^)  which  correspond  to  those 
across  a normal  shocJc  wave  (see  Appendix  A) . 
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In  order  to  establish  the  validity  of  solutions  and  to  evaluate  the 
pressure  ratio  (P2/P2)  for  a non-uniform  flow,  it  is  necessary  to  consider 
the  entropy  equation.  It  is  convenient  to  consider  the  entropy  at  Stations 
1,  2 and  3 using  the  free  stream  entropy  as  a base.  Therefore 


S,  - S = m (s  , - s)  + m.  (s, 
1 “ P pi  “ 1 1 


s_) 
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and  at  Station  2 
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dy 
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which  upon  integration  becomes,  (see  Appendix  B for  details) 
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The  Second  Law  of  Thermodynamics  applied  to  the  mixing  process  dictates 


that  S^  > 
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The  pressure  ratio  P2/P3  <=  'P2^'Pa)  ^ expressed  in  terras  of  the  entropy 

change  between  Stations  2 and  3,  using  Equations  43  and  44,  as  follows. 


--(-rlr  \r'') 


Since  = T^2'  mass  flow  averaged  enthalpy  relationship  described  by 
Equation  27  can  be  expressed  at  Station  3 in  the  form. 


1 + 3 — +2  - 
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dividing  Eq.  27  by  Eq.  46 


p,  V p2  S'  L ^ s P2  S'- 


Combination  of  Equations  45  and  47  to  eliminate  P2/P3  rearrangement  to 
solve  for  A^,  yields  the  result 
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where  P2/P00  evaluated  from  Equation  28. 


Elimination  of  P2/P00  fro«®  Equation  45  and  47  provides  the  expression 
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which  is  in  a convenient  form  for  use  in  evaluating  the  required  effective  outlet 
area  ratio  (6*)  since  from  continuity  applied  between  Stations  2 and  3,  similarly 
to  Equation  20, 
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or,  since  from  Equation  48 
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SECTION  III 

PERFORMANCE  PARAMETERS 


1.  Thrust  Augmentation  ((t>) 

Thrust  augmentation  is  defined  as  the  ratio  of  the  net  thrust  of  an  ejector 
to  the  net  thrust  of  a reference  "free"  jet  whose 


a)  Jet  power  is  equal  to  that  of  the  energized  flow  of  the  ejector 

b)  Mass  flow  is  equal  to  that  of  the  energized  flow  of  the  ejector 

c)  Ratio  of  mechanical  energy  to  thermal  energy  is  equal  to  that  of 
the  energized  mass  flow  of  the  ejector. 


Under  these  constraints,  the  net  thrust  of  the  reference  jet  is 
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2. 


Propulsive  Efficiency  ( n ) 


The  propulsive  efficiency  of  a free  jet  is  commonly  described  as  the  ratio 
of  the  useful  energy  to  the  mechanical  energy.  Thus  in  general 


Thrust  X Velocity 
Mechanical  Energy 


In  the  case  of  an  ejector 
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(61) 


since  the  mechanical  energy  is  that  of  its  reference  free  jet  whose  mass  flow, 
jet  power  and  ratio  of  mechanical  to  thermal  energy  are  equal  to  those  of  the 
ejector's  energized  flows,  as  in  the  reference  jet  used  to  describe  thrust  augmen- 
tation. 

For  non-air-breathing  engines  (roc)cets)  , 


n, 


p , rocl^et 


1 + (uyv.) 


q-,  and,  n . = 4)  n , ^ 
2 P,eD  p,roc)cet 


(62) 


3.  Relative  Power  (R  ) 

E_ 


It  is  frequently  of  interest  to  compare  thrusters  on  the  basis  of  their  power 
requirements  when  their  net  thrusts  are  equal.  Since  the  equality  of  net  thrust 
does  not  uniquely  define  the  power  relationship  among  thrusting  systems  it  is 
necessary  to  impose  additional  constraints.  In  a system  utilizing  compressible 
fluids  two  additional  constraints  are  required,  and  for  purposes  of  the  comparison 
of  ejector  systems  with  free  jets  these  constraints  are  chosen  to  be  the  equality 
of  energized  mass  flow  and  the  equality  of  the  ratio  of  mechanical  to  thermal 
energy.  Thus 


Ejector  Power 
Reference  Jet  Power 
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which  in  view  of  the  constraints  described  above,  reduces  to  the  form 
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The  equality  of  net  thrusts  and  energized  mass  flows  can  be  expressed  as 


i{)(V  - U ) = V ^ - U 

ptu  to  poo , r e f 


or 


V + u = (J)(v  + U ) - 2(<j)  - 1)U 


and  therefore 


(V  + u ) 

p«’  ■” 


(65) 


+ UJ  - - DU 


- (4)  - i)n. 


- (})(())  - Dn, 


P»e3  p.ref 

as  indicated  by  a comparison  of  Equations  60  and  61, 


since  h . = tjiTi 

p,ej  p.ref 

As  can  be  observed  from  Equation  65,  the  ejector  can  be  utilized  to  produce 


a given  net  thrust  with  a considerably  smaller  required  jet  power  than  that 
required  by  a free  jet  from  a gas  generator  whose  net  thrust,  mass  flow  and  ratio 
of  mechanical  to  thermal  energies  are  equal  to  those  of  the  ejector's  injected 


jets.  For  example  when  V = U , R = l/4>  and  as  in  most  application  V >>  U^  and 
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4 . Efficiencies 

To  describe  the  performance  of  combination  thrusters  in  which  various  types 
of  gas  generators  are  used  in  conjunction  with  ejectors,  it  is  necessary  to 
evaluate  the  influence  of  the  ejector  upon  the  cycle  efficiencies  of  the  gas 
generator.  Thus  conventionally,  the  propulsive  efficiency  of  a free  jet  is 


B Useful  Energy 
^ Mechanical  Energy 
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and  for  an  ejector,  since  the  propulsive  efficiency  referred  to  the 

injected  mecheuiical  energy 
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since  the  useful  energy  of  the  combination  thruster  is  (p  times  the  useful  energy  of 

its  reference  "free"  jet  which  has  the  identical  mechcmical  energy  as  that  of  the 

energized  jet  of  the  ejector.  This  definition  of  propulsive  efficiency  for  ejectors 

is  more  useful  than  that  which  describes  the  ejector  as  a free  jet.  As  a free 

jet  if  the  ejector  exhaust  has  a velocity  U , its  propulsive  efficiency  is 
2 U 
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P.ej 


(68) 


since  mechanical  energy  is  not  conserved  in  the  ejector  process. 

Evaluation  of  propulsion  system  performance  by  means  of  the  propulsive 
efficiency  is  unsatisfactory  for  comparison  of  various  air  breathing  gas  gener- 
ator/ejector combination  systems  since  defined  in  this  mcuiner  useful  energy  is 
compared  only  to  the  jet  mechanical  energy;  the  thermal  energy  is  neglected. 
Therefore  the  overall  efficiency  (n^)  is  more  appropriate  for  such  comparisons. 
The  overall  efficiency  is  defined  as  the  ratio  of  useful  energy  to  the  total  jet 
energy 
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Useful  Energy 
Total  Jet  Energy 
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for  a "free"  jet. 

As  in  the  definition  of  propulsive  efficiency,  the  overall  efficiency  of  a 
combination  gas  generator/ejector  system  is 
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since  the  useful  energy  is  ())  times  that  of  the  free  reference  jet  and  the  total 
jet  energy  of  the  ejector  is  equal  to  that  of  its  free  reference  jet. 
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Since  the  thermal  efficiency  (H^)  of  a system  is  defined  conventionally  as 
the  ratio  of  mechanical  energy  to  total  energy,  it  can  be  expressed  for  a "free 
jet  as 


I _ Mechanical  Energ; 
t Total  Jet  Energy 


— (71) 

P P op, ref  o'” 

This  definition  of  thermal  efficiency  represents  the  distribution  of  the 
jet  energy  content,  rather  than  the  thermodynamic  cycle  efficiency.  For 
example,  an  Isentropically  compressed  jet  (fzmjet)  contains  100%  mechanical 
energy  (100%  conversion  of  mechanical  energy  to  jet  energy  from  any  power 
source,  such  as  jet  engines,  reciprocating  engines,  or  electric  motors),  and 
has  a 100%  thermal  efficiency  as  defined  by  Equation  71. 

For  an  ejector  in  which  the  mass  flow  is  augmented,  the  thermal  efficiency  is 


T)  . may  be  larger  than  or  smaller  than  n depending  upon  the  magnitude  of  the 
t , e]  t 

term  + Uoo^/^^poo  • Examples  illustrating  the  magnitude  of  the  ratio 

(mechanical  energy  output  of  ejector /mechanical  energy  output  of  refer 

t , ej  t 

ence  jet)  are  presented  in  the  following  section. 

Since  for  a free  jet. 


by  definition,  the  overall  efficiency  of  a combined  gas  generator/ejector  system  is 


o,ej  o 

indicating  that  the  gas  generator/ejector  combination  can  produce  times  the  useful 

work  of  that  which  can  be  produced  by  its  reference  "free"  jet  by  virtue  of  its 

ability  to  convert  lost  mechanical  and  thermal  energy  to  useful  energy. 

It  should  be  emphasized  that  n . since  the 

o,ej  'p,ej  't,e:  ^ 'p,ej  't,ej 

mechanical  energy  input  does  not  equal  the  mechanical  energy  output  of  the 

ejector  due  to  its  intrinsic  change  of  mechanical  energy  during  mixing  and  the 

definition  of  ri  . discussed  previously. 

P»e: 


k 


SECTION  IV 


EJECTOR  PERFORMANCE 


The  mixing  of  primary  and  induced  flows  in  an  ejector  can  result  in  two 
distinctly  different  conditions  at  the  conclusion  of  the  mixing  process. 
Mathematically  this  is  indicated  by  the  quadratic  nature  of  the  equation  for 
(Equation  34) . Physically  this  dual  solution  represents  two  different  ejector 
geometries.  It  is  shown  in  Appendix  A that  after  complete  mixing  the  Mach  No.  of 
the  flow  is  less  than  1.0  when  using  the  negative  sign  in  the  solution  of  the 
quadratic  equation  and  greater  than  1.0  when  the  positive  sign  is  used,  emd  it 
is  proven  that,  in  the  absence  of  skin  friction,  the  two  Mach  No's  are  related 
by  a normal  shock  wave  relationship  at  the  end  of  complete  mixing. 

The  first  solution  corresponding  to  M2  <1  always  represents  a flow  in 
which  the  total  entropy  change  is  greater  tham  zero  (S2  m^  > s^^  m^  + s^^^  m^) 
and  which  therefore  is  consistent  with  the  Second  Law  of  Thermodynamics.  The 
second  solution  is  always  consistent  with  the  Second  Law  of  Thermodynamics  when 
Mj^  > 1 (supersonic  mixing),  but  at  low  subsonic  values  of  Mj^,  it  can  result  in 
mathematically  correct  solutions  which  represent  a decrease  of  entropy  or  in  which 


(s_  m < s,  m.  + s , m ) across  the  mixing  process,  emd  which  therefore  is  in 
2cliplp 

violation  of  the  Second  Law  of  Thermodynamics.  It  therefore  is  essential  to 
evaluate  the  entropy  change  before  utilizing  the  results  of  these  calculations. 

A further  characteristic  of  the  solutions  of  the  equations  for  the  flow  of  a 
compressible  fluid  through  an  ejector  is  evidenced  by  regions  where  the  solutions 
are  complex.  The  limits  of  these  regions  are  those  points  at  which  the  discriminant 
in  the  solution  of  Equation  34  is  zero.  These  t\ra  values  represent  a subsonic 
and  a supersonic  induced  flow  at  Station  1.  Values  of  Mj^  between  these  limits 
produce  complex  solutions  which  indicate  that  the  flow  has  become  choked  (M2  ° 1.0) 
at  or  prior  to  the  completion  of  mixing.  At  the  choking  limits,  the  flow  parameters 
are  Identical  for  each  of  the  two  solutions  and  the  precise  values  of  M^  at  which 
they  occur  depend  upon  the  properties  of  the  Injected  gas,  the  ejector  geometry 
and  the  flight  Mach  number. 

The  solutions  to  these  equations  are  also  useful  in  the  selection  of  ejector 
flow  characteristics  which  produce  optimum  performance.  In  the  regions  where 
valid  solutions  are  obtained,  there  are  in  general  three  distinct  points  or  flow 
conditions  which  represent  the  best  performance. 
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Extensive  numerical  analysis  over  large  ranges  of  flight  Mach  No's,  and 
injected  gas  characteristics  have  disclosed  that  in  general  the  first  solution 
represents  a flow  in  which  the  meucimum  thrust  augmentation  occurs  at  subsonic 
values  of  The  second  solution  represents  a flow  in  which  the  maximum  thrust 

augmentation  exists  at  values  of  greater  than,  or  slightly  less  than  1.0,  and 
a limiting  value  of  thrust  augmentation  at  a subsonic  value  of  Mj^.  This  limiting 
value  of  thrust  augmentation  occurs  at  a point  where  the  total  entropy  of  the 
mixed  flow  at  Station  2 is  equal  to  the  sum  of  the  entropies  of  the  unmixed 
primary  and  induced  flows. 

Further  characterization  of  the  functional  dependence,  limitations  auid 

optimization  of  the  thrust  augmentation  under  both  solutions  are  discussed  in 

terms  of  two  specific  examples  to  provide  eui  indication  of  the  method  for  optimal 

design  of  translating  ejectors  in  a compressible  medium.  These  calculations  assume 

that  the  mixing  is  completed  at  Station  2,  and  that  the  flow  is  isentropic  between 

Stations  2 and  3 (s,“  3_) , emd  that  n„  = 1,  C, .=  C = 0,  6*  = 6,  U = u , y = 1.4. 

j z N di  F e 00 

In  addition,  the  calculations  were  made  for  a fixed  value  of  a to  avoid 

00 

variations  of  ejector  size  compared  to  primary  jet  size  with  variations  of  M^. 

1.  Low  speed  example 

During  low  speed  flight,  combination  gas  generator/ejector  thrusters 
display  good  overall  system  characteristics  when  the  gas  generator  is  of  the 
fanjet  type  or  when  the  ejector  is  energized  by  a fluid  which  is  mechanically 
compressed  (isentroplcally)  without  the  addition  of  heat  energy.  This  is  a 
result  of  the  fact  that  at  low  flight  speeds,  high  temperature  jets  have  low 
thermal  efficiencies.  An  example  of  the  characteristics  of  the  curves  representing 
the  ejector  performance,  flow  properties,  ejector  geometry  and  thermodynamic  cycles 
for  both  solutions  during  subsonic  flight  speeds  are  illustrated  on  Figure  2 . 

a.  First  solution 

As  illustrated  on  Figure  2,  the  first  solution  has  a maximum  thrust 
augmentation  at  a value  of  close  to  1.0,  at  which  point  the  transfer  of 
mechanical  energy  from  the  reference  jet  to  the  ejector  flow  is  also  a maximum. 

The  flow  characteristics  illustrated  on  Figure  2 indicate  a maximum  entrainment 
ratio  and  a minimum  (but  positive)  entropy  change  during  mixing.  Temperatures 
and  pressures  at  various  stations  within  the  ejector  are  also  presented  in  a 
dimensionless  form  on  this  figure.  The  local  Mach  No's  and  the  corresponding 
area  ratios  are  also  illustrated  on  Figure  2 , where  it  can  be  obs  srved  that  under 
the  first  solution  the  value  of  M^^  at  the  optimal  point  is  large  than  the  free 


stream  Mach  No.,  thus  requiring  an  accelerating  inlet.  The  fact  that  is 
smaller  than  and  that  both  are  subsonic  indicates  a decelerating  outlet  or 
diffuser  is  required  at  the  specified  conditions,  resulting  in  a conventional 
ejector  design. 

The  Mollier  Diagram  illustrated  on  Figure  2 at  this  optimal  point  of  the 
first  solution  describes  the  properties  of  the  energized  flow  and  the  ejector 
mixing  process.  Since  the  primary  injected  gas  emd  the  induced  flow  are  assumed 
to  undergo  isentropic  processes  from  their  undisturbed  state  to  Station  1,  in 
this  example,  their  specific  entropies  at  Station  1 are  equal,  thus  s^^  = s^^^  “ ®oo 
euid  T^  = Mixing  in  this  case  results  in  aui  increased  entropy  for  both  streams 

and  is  concluded  with  a small  increase  in  specific  enthalpy  at  the  point,  where 
primary  and  induced  flows  are  completely  mixed.  The  recompression  to  ambient 
pressure  is  then  assumed  to  be  isentropic  as  indicated. 

b.  Second  solution  (supersonic  mixing) 

The  second  solution  presents  two  distinct  flow  configurations  which  represent 
performance  characteristics  of  interest.  One  of  these  occurs  with  supersonic  mixing 
(Mj^  > 1)  and  as  shown  on  Figure  2,  represents  a situation  in  which  the  overall 
mechanical  energy  of  the  ejector  flow,  is  a larger  percentage  of  the  mechauiical 
energy  of  the  reference  jet  than  that  at  the  optimal  point  of  the  first  solution. 

The  entropy  change  during  mixing,  the  entrainment  ratio  and  the  pressures  and 
temperature  within  the  ejector  are  also  described  on  Figure  2 for  comparison  with 
the  first  solution.  As  indicated,  under  the  second  solution,  this  optimal  point 
has  a positive  but  smaller  entropy  change  and  a smaller  entrainment  ratio  than 
that  of  the  optimal  condition  under  the  first  solution.  The  Mach  No's,  and  corre- 
sponding area  ratios  for  this  condition  are  presented  on  Figure  2 where  it  is 
evident  that  the  inlet  is  accelerating  and  the  outlet  is  decelerating.  However 
since  is  subsonic  and  Mj^  is  supersonic,  the  ejector  inlet  must  have  a converging/ 
diverging  De  Lavel  Nozzle,  and  since  M2  is  supersonic,  and  M^  is  subsonic,  the 
ejector  outlet  must  be  a converging/diverging  De  Laval  Nozzle.  The  Mollier  Diagram 
under  this  condition  is  presented  on  Figure  2. 

Some  additional  features  associated  with  this  supersonic  mixing  point  under 
the  second  solution  include  the  facts  that  the  thrust  augmentation  is  only  slightly 
dependent  upon  or  outlet  area  ratio  (6)  in  this  region  and  that  the  improved 
performance  occurs  with  a smaller  are&  ratio  (6)  than  that  of  the  first  solution. 


r'  r 


c.  Second  solution  (limit  of  subsonic  mixing) 

The  other  point  of  interest  under  the  second  solution  provides  a limiting 
performance  within  the  Second  Law  of  Thennodyn2unics  which  occurs  during  subsonic 


mixing  < 1) . As  indicated  on  Figure  2,  this  point  represents  the  largest 


thrust  augmentation  among  all  three  configurations  of  interest,  a transfer  of 
100%  of  the  mechanical  energy  of  the  reference  jet  to  the  ejector  and  as  shown 
on  Figure  2,  zero  total  entropy  change  during  mixing.  The  entrainment  ratio 
is  not  a maiximum  at  this  point  but  is  very  close  to  that  of  the  first  solution. 
Temperatures  and  pressures  are  also  indicated  on  Figure  2.  The  inlet  of  the 
ejector  at  this  condition  is  accelerating  and  the  outlet  is  a converging/diverging 
De  Laval  Nozzle  providing  a mecins  for  reduction  of  the  Mach  No.  from  a supersonic 
value  at  Station  2 to  a subsonic  value  at  Station  3 as  indicated  on  Figure  2. 

The  Mollier  Diagram  for  this  cycle  is  presented  on  Figure  2 where  it  is 
shown  to  degenerate  to  a single  isentropic  line  as  a result  of  the  assumption 
of  isentropic  compression  and  expansion  and  the  selection  of  the  point  where  the 
total  entropy  of  the  mixed  flow  is  equal  to  the  sum  of  the  total  entropies  of  the 
primary  and  induced  flows.  The  loss  of  mechanical  energy  within  the  ejector  in 
this  case  is  zero. 

It  is  of  interest  to  note  (Figure  2)  that  at  this  limiting  point  of  subsonic 
mixing  under  the  second  solution  the  propulsive  efficiency  of  the  combined  fanjet/ 
ejector  system  can  achieve  values  in  the  vicinity  of  82%  compared  to  44%  for  the 
fanjet  alone.  This  indicates  a recovery,  by  the  ejector,  of  a large  percentage  of 
the  kinetic  energy  which  is  lost  by  the  fcmjet  when  operating  as  a free  jet. 
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2.  High  Speed  Example 

At  high  flight  speeds,  the  equations  describing  the  laws  of  mass  flow  and 
energy  conservation  and  the  momentum  theorem  produce  two  distinct  solutioi s as  in 
the  case  of  low  speed  flight.  However,  at  high  speeds,  the  use  of  high  temperature 
injected  gas  becomes  practical  from  the  point  of  view  of  gas  generator  thermal 
efficiency  and  this  leads  to  a limitation  of  ejector  performance  due  to  thermal 
choking.  This  phenomenon  occurs  when  the  Mach  No.  of  the  completely  or  partially 
mixed  flow  reaches  a value  of  1.0,  and  requires  the  avoidance  of  a set  of  design 
configurations  corresponding  to  values  of  near  1.0.  This  region  is  depicted  on 
Figure  3,  as  a reinge  of  values  of  bounded  by  points  described  as  the  choking 
limit.  Other  characteristics  of  the  solutions  to  the  flow  equations  at  high  speeds 
are  described  in  terms  of  the  performance  eind  flow  properties  of  a ramjet/ejector 
in  regions  outside  the  choking  limit,  at  a flight  Mach  No.  of  2.0. 

a.  First  Solution 

As  before,  the  first  solution  is  that  solution  to  the  equations  which  result 
in  a subsonic  flow  after  complete  mixing.  The  aerodynamic,  thermodynamic  and 
geometric  characteristics  of  the  ejector  at  supersonic  flight  speeds  are  described 
on  Figure  3. 

As  illustrated  on  Figure  3,  the  thrust  augmentation  with  subsonic  mixing 
and  subsonic  mixed  flow  after  complete  mixing  has  a maximum  value  at  which  the 
transfer  of  mechanical  energy  from  the  reference  jet  to  the  ejector  is  99%,  and 
the  propulsive  efficiency  for  this  case  achieves  a value  of  96%,  with  an  overall 
efficiency  of  43%  while  under  the  same  conditions,  the  ramjet  itself  achieves  a 
propulsive  efficiency  of  56%  with  an  overall  efficiency  of  only  25%.  Therefore 
it  is  apparent  that  a portion  of  the  waste  energy  of  the  ramjet  is  converted  to 
useful  energy  by  the  ejector. 

At  this  optimal  point  with  subsonic  flow  at  the  start  and  at  the  conclusion 
of  complete  mixing,  the  total  entropy  change  is  positive,  signifying  compliance 
with  the  Second  Law  of  Thermodynamics  as  illustrated  on  Figure  3.  Further  as 
shown  on  this  figure,  the  entrainment  ratio  is  not  at  its  maximum  nor  is  the 
entropy  production  at  its  minimum.  Pressure  euid  temperature  at  the  various 
stations  within  the  ejector  are  also  described  on  this  figure. 


t 
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The  Mach  No.  and  area  ratios  of  the  ejector  flows  are  described  on  Figure  3. 

As  illustrated,  operation  at  this  optimal  point  requires  a converging/diverging, 
decelerating  De  Laval  Nozzle  at  the  inlet  to  the  ejector  and  a converging/diverging, 
accelerating  supersonic  nozzle  at  the  exit  of  the  ejector.  The  area  ratio  of  the 
converging/diverging,  accelerating  exit  nozzle  is  0.94,  somewhat  less  than  1.0. 

The  Mollier  Diagram  showing  the  specific  enthalpy  vs.  specific  entropy  changes 
occuring  in  the  primary,  energized  flow  and  in  the  induced  and  mixed  flows  are 
presented  on  Figure  3.  The  large  decrease  in  specific  entropy  and  specific 
enthalpy  of  the  primary  flow  and  the  corresponding  small  increase  of  these 
properties  in  the  induced  flow  result  in  a net  increase  of  the  total  entropy 
after  mixing  as  a consequence  of  the  large  entrainment  ratio.  Thus  the  cycle  at 
this  optimal  point  is  not  in  violation  of  the  Second  Law  of  Thermodynamics. 

b.  Second  Solution  (supersonic  mixing) 

The  second  solution,  as  before  at  low  speeds,  displays  two  distinct  flow 
configurations  of  interest.  One  of  these  occurs  with  supersonic  mixing  and  in  this 
particular  example  is  coincident  with  the  first  solution  since  the  maximum  thrust 
augmentation  lies  exactly  at  the  supersonic  choking  limit. 

As  illustrated  on  Figure  3,  the  thrust  augmentation  and  the  transfer  of 
mechanical  energy  are  smaller  than  those  with  subsonic  mixing  but  the  advantage  may 
result  from  the  possibility  of  reduced  inlet  loss  with  supersonic  values  of 
Propulsive  and  overall  efficiencies  remain  high  compared  to  those  of  the  rarr.jet 
alone. 

Figure  3 indicates  a positive  change  of  total  entropy  and  an  entrainment 
ratio  close  to  the  m«u{imum  for  this  condition.  Pressures  and  temperatures  through- 
out the  ejector  are  also  described  on  this  figure. 

The  local  Mach  No's,  euid  area  ratios  are  described  on  Figure  3 where  as 
shown  the  inlet  is  a converging,  decelerating  nozzle  (Mj^  < M^,  but  both  are  super- 
sonic) , while  the  outlet  is  a diverging,  accelerating  nozzle.  The  overall 
appearance  of  this  ejector  is  identical  to  low  speed,  low  temperature  ejectors  of 
the  first  solution. 

The  Mollier  Diagram  shown  on  Figure  3,  is  simlar  to  that  with  subsonic 
mixing,  but  as  a consequence  of  the  larger  entrainment  ratio,  the  primary  fluid 
has  a larger  decrease  of  specific  entropy  while  the  entrained  flow  has  a smaller 
increase  of  specific  entropy. 
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c.  Second  Solution  (limit  of  subsonic  mixing) 


This  point  is  a limit  imposed  by  the  Second  Law  of  Thermodynamics  at  which 
the  change  of  total  entropy  is  zero,  and  where  the  thrust  augmentation  and  the 
transfer  of  mechanical  energy  are  considerably  larger  than  at  the  other  optimal 
conditions . 

As  illustrated  on  Figure  3,  the  large  value  of  thrust  augmentation  at  the 
point  where  the  performance  is  limited  by  the  Second  Law  is  a result  of  the  con- 
version of  thermal  energy  from  the  energized  flow  to  mechanical  energy  of  the 
ejector  flow.  Propulsive  efficiencies  in  excess  of  100%  are  achieved  under  these 
conditions  while  overall  efficiencies  are  extremely  large  but  remain  less  than 
100%.  Thus  there  is  no  violation  of  the  law  of  energy  conservation. 

As  indicated  on  Figure  3,  the  total  entropy  change  is  zero  and  the 
entrainment  ratio  is  not  at  its  maximum  under  these  conditions.  The  pressure  and 
temperature  at  Station  2 are  at  a minimum  at  this  condition  and  as  shown  on 
Figure  3,  the  Mach  No.  at  Station  2 is  a maximum.  Thus  the  high  Mach  No.  at 
Station  2 is  achieved  primarily  by  a static  pressure  reduction  (or  expansion) . 

Figure  3 illustrates  the  Mach  No's,  and  area  ratios  required  to  achieve  the 
flow  characteristics  previously  described.  As  indicated  the  inlet  must  be  of  the 
converging/diverging,  decelerating,  De  Laval  Nozzle  with  subsonic  induced  flow  at 
Station  1.  The  outlet  which  accelerates  the  flow  from  a supersonic  value  at 
Station  2 to  a higher  supersonic  value  at  Station  3 is  diverging  cind  accelerating. 

The  Mollier  Diagram  for  operation  at  this  point  is  illustrated  on  Figure  3. 

It  is  somewhat  similar  to  the  corresponding  diagram  for  the  other  optimal  points 
at  this  conditions  but  since  the  change  of  total  entropy  is  smaller  (zero)  the 
decrease  of  specific  entropy  of  the  energized  fluid  is  larger  and  the  increase  of 
specific  entropy  of  the  induced  flow  is  smaller  in  magnitude  than  those  of  the 
other  optimal  points,  despite  an  entrainment  ratio  somewhat  smaller  than  that 
corresponding  to  supersonic  mixing. 
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Figure  3.  Ejector  in  High  Speed  Flight 
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Propulsive  Efficiency 


en^jerature,  T/T  Temperature,  T/T  Temperature,  T/T 


It  Is  interesting  to  note  that  ejectors  operating  underwater  or  zmalyzed 
with  the  use  of  incompressible  flow  theory  (References  2 amd  3)  always  require 
accelerating  (convergent)  inlets  emd  decelerating  (divergent  or  diffusing)  outlets 
to  achieve  optimal  performance  under  the  single  realistic  solution  which  exists. 

In  contrast,  the  introduction  of  compressibility  effects  (variable  temperature 
and  density)  into  the  analysis  results  in  the  existence  of  two  distinct,  realistic 
solutions  representing  ejectors  which  may  require  decelerating  inlets  and 
accelerating  outlets.  Depending  upon  the  flight  Mach  No.  emd  the  characteristics 
of  the  injected  gas,  the  ejector  will  require  one  of  nine  configurations.  These 
nine  configurations  are  the  result  of  the  possibility  of  a choice  among  three 
types  of  inlets  (converging,  diverging  or  De  Laval  Nozzle)  and  three  types  of 
outlets  (converging,  diverging  or  De  Laval  Nozzle) . 

The  first  of  the  two  solutions  to  the  compressible  flow  equations  results 
in  ejector  conflgiiratlons  2uid  flow  patterns  which  at  low  flight  speeds,  powered 
by  low  temperature  jets  closely  resemble  those  which  are  derived  from  the  theory 
of  incompressible  flow.  At  high  flight  speeds  or  high  temperature  primary  jets 
this  solution  indicates  a radical  departure  from  the  ejector  configurations  derived 
from  incompressible  flow  theory.  High  thrust  augmentation  can  be  achieved  under 
this  first  solution  over  the  entire  range  of  flight  speeds  provided  the  ejector 
geometry  emd  properties  of  the  injected  gas  are  properly  selected,  as  will  be 
shown  on  Figures  6 emd  7. 

The  second  solution  with  supersonic  mixing  and  supersonic  flight  speeds 
requires  an  ejector  configuration  which  geometrically  resembles  those  derived 
from  incompressible  flow  theory,  but  which  perform  the  opposite  function  due  to 
the  difference  between  the  supersonic  eind  subsonic  area/Mach  No.  relationship. 

At  low  flight  speeds  the  geometry  of  the  ejector  differs  considerably  from 
those  derived  under  incompressible  flow  theory,  but  the  function  can  be  similar  to 
those  derived  under  incompressible  flow  theory  if  powered  by  low  temperature  jets. 
With  proper  selection  of  the  ejector  geometry  and  the  properties  of  the  injected 
gas,  the  ejectors  developed  with  the  use  of  the  second  solution  can  produce  better 
ejector  performance  them  those  derived  under  the  first  solution. 

The  importance  of  a careful  selection  of  ejector  geometry  and  injected  gas 
characteristics  cannot  be  overemphasized  if  high  ejector  performance  is  to  be 
achieved  during  translatory  conditions,  under  either  of  the  two  solutions  to  the 
flow  equations. 


The  general  features  of  these  tvfo  solutions  have  been  described  in  terms 
of  two  examples  at  two  discrete  operational  conditions.  It  has  been  shown  that 
In  general,  for  a given  condition,  there  exist  three  types  of  ejector  or  flow 
configurations  which  are  of  Interest  In  the  design  of  ejectors.  Two  of  these 
configurations  are  the  result  of  analysis  using  the  second  solution.  One  of 
these  two  exists  when  the  induced  flow  is  subsonic  at  the  start  of  mixing,  while 
the  other  exists  when  the  Induced  flow  Is  supersonic  (or.  In  some  special  cases, 
slightly  less  than  sonic)  at  the  start  of  mixing.  The  third  configuration  of 
interest  Is  the  result  of  the  analysis  when  the  first  solution  Is  used  and  the 
induced  flow  Is  subsonic  at  the  start  of  mixing. 

Mollier  diagrams  illustrating  the  specific  enthalpy  as  a function  of  the 
specific  entropy  for  each  of  the  selected  design  configurations  are  presented 
to  describe  the  processes  undergone  by  the  primary,  induced  and  mixed  flows 
in  the  ejector.  These  diagrams  are  particularly  useful  in  describing  modifi- 
cations to  the  basic  gas  generator  thermodynamic  cycle  resulting  from  the  use 
of  ejector  thrust  augmentation. 

To  illustrate  the  feasibility  of  ejector  propulsion,  the  thrust  augmen- 
tation and  propulsive  efficiency  of  the  reference  jet  have  been  shown  for 
representative  gas  generator/ejector  combination  thrusters  operating  at  the 
three  design  points  over  a range  of  flight  Mach  No's  in  the  following  section. 


3.  Combination  Thrusters 


The  performance  of  combination  thrusters  over  a range  of  flight  Mach  No's. 
fr<»n  0 to  3 has  been  evaluated  for  fan jet,  ramjet,  turbojet  and  rocket  type  gas 
generators.  Each  type  of  gas  generator  is  of  interest  for  different  applications 
and  under  different  operational  conditions,  depending  upon  the  requirements  of 
size,  efficiency,  etc.  For  example,  the  ideal  performance  of  an  ejector  generally 
improves  slowly  with  increasing  values  of  the  ratio  of  ejector  mixing  duct  size 
(A^)  to  primary  jet  (expanded  to  ambient  pressure)  size.  In  the  examples 

described  below,  this  ratio  (a  ) is  maintained  at  a fixed  value  of  20,  for 

00 

discernment  of  the  trend  of  performance  as  a function  of  flight  Mach  No.  for 
each  type  of  gas  generator/ejector  combination. 

The  calculation  of  performance  of  combination  ejector  thrusters  is  based 
upon  fixed  values  of  AP/p^,  the  ratio  of  pressure  rise  (P^^.  “ to  the  ambient 

pressure,  euid  AT/T  , the  ratio  of  temperature  rise  (T  - T ) to  the  ambient 

oo  * Op  0^0 

temperature.  In  this  manner  the  stagnation  properties  of  the  primary  jet  remain 
realistic  with  respect  to  flight  Mach  No.  eind  altitude. 


The  calculation  of  the  fan jet/e jector  performance  is  based  upon  the 
assumption  that  the  compression  is  isentropic  as  previously  described.  Thus 
since  the  energized  fluid  contains  only  mechanical  energy  the  overall  efficiency 
is  equal  to  the  propulsive  efficiency  or  in  other  words  the  thermal  efficiency 
as  defined  in  this  document  is  100%. 

Figure  4 presents  the  thrust  augmentation  over  a range  of  flight  Mach  No's, 
from  zero  to  three  when  AP/p^  = 10  and  “ 20.  This  represents  one  particular 
condition  at  which  this  device  can  operate  and  as  indicated  the  thrust  augmentation 
decreases  rapidly  with  increasing  flight  Mach  No.  with  all  three  types  of  optimal 
ejector  configurations.  However  the  ejector  designed  with  subsonic  mixing  under 
the  second  solution  provides  superior  performance  to  the  other  configurations 
over  the  entire  range  of  operational  conditions.  As  indicated,  it  is  essential  to 
use  the  design  dictated  by  the  second  solution  as  a design  criterion  for  transonic 
or  supersonic  ejectors  powered  by  fanjets.  Successful  application  of  a fanjet 
ejector  requires  careful  selection  of  the  fanjet  pressure.  For  example,  at  flight 
Mach  Number  of  0.3,  the  best  remge  of  AP/p^  for  all  three  characteristic  design 
points  lies  between  1 and  4 as  shown  on  Figure  5 for  an  inlet  area  ratio  (a^) 
of  20. 


b.  Ramjet/Ejectors 


The  ramj at/ejector  combination  thruster  is  a system  in  which  the  injected 
gas  (efflux  from  a rtunjet)  contains  large  quantities  of  thermal  emd  mechanical 
energy  produced  as  a result  of  the  Brayton  Cycle  with  compression  due  to  the 
motion  of  the  system,  with  Ap/p^  = 0,  and  under  the  other  stated  conditions, 
the  performance  of  the  systmn  is  presented  on  Figure  6.  As  indicated,  with 
subsonic  mixing  this  system  provides  extremely  high  ejector  performance  over  the 
entire  range  of  flight  Mach  No's.,  for  ejector  designs  derived  under  both  solutions. 
Supersonic  mixing  provides  high  performeuice  only  at  very  large  flight  Mach  No's. 

Since  the  thermal  efficiency  is  small  at  low  flight  speeds  the  overall 
efficiency  is  also  small  and  therefore  this  combination  is  of  little  interest  as 
a low  cruise  speed  propulsive  system.  However  the  simplicity  cind  reliability  of 
the  ramjet  as  a gas  generator  provides  cm  attractive  system  for  high  speed 
propulsion  auid  in  combination  with  an  ejector,  this  system  can  be  utilized 
effectively  during  the  accelerating  low  speed  portion  of  the  flight  spectrum 
of  a high  cruising  speed  vehicle. 

At  high  flight  speeds,  the  high  ejector  thrust  augmentation  shown  on 
Figure  6 for  At/T^  «•  10,  depends  strongly  upon  the  temperature  and  the  design 
configuration.  For  example  using  the  first  solution  design  point  with  subsonic 
mixing  it  is  shown  on  Figure  7 that  the  thrust  augmentation  depends  strongly 
upon  AT/T^.  With  decreasing  values  of  At/T^  the  thrust  augmentation  drops  off 
rapidly  over  the  entire  range  of  flight  speeds.  Using  the  second  solution  design 
point  with  supersonic  mixing  results  in  a similar  dependence  of  thrust  augmentation 
upon  the  temperature  (Figure  7) , but  good  performance  is  achieved  only  at  high 
flight  speeds.  This  limitation  at  low  speeds  for  the  second  solution  with 
supersonic  mixing  is  of  little  consequence  however  since  as  indicated  above,  the 
ramjet/ejector  is  not  a practical  system  for  cruise  at  these  low  speeds.  With 
subsonic  mixing  the  second  solution  also  requires  high  temperature  primary  jets 
for  high  ejector  performance.  Thus  it  is  evident  that  regardless  of  the  design 
point  utilized,  the  ramjet/ejector  system  operates  effectively  at  high  speeds 
and  high  temperatures.  If  low  speed  flight  with  this  power  plant  is  required, 
effective  performance  can  only  be  achieved  with  subsonic  mixing. 
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c.  Turbojet/Ejector 

Figure  8 Illustrates  the  performance  of  a turbojet/ejector  combination 
system  over  a remge  of  flight  Mach  No's,  from  0 to  3 with  injected  gas  having 
high  mechanical  and  thermal  energies  (AP/p^  • 10,  AT/T^  = 10) . As  indicated,  at 
high  flight  speeds  the  thrust  augmentation  is  improved  over  the  fanjet/e jector 
and  is  equivalent  to  the  ramjet/ejector  under  these  conditions.  At  low  speeds 
the  turbojet/ejector  having  the  stated  properties  of  the  injected  gas  is  competi- 
tive with  the  fan jet  under  all  three  design  conditions. 

At  intermediate  flight  Mach  Numbers  (near  0.5),  this  turbojet/ejector  having 
the  properties  stated  above  produces  very  low  ejector  performance,  if  designed 
according  to  the  first  solution,  due  to  the  fact  that  the  design  requirements  for 
the  use  of  hot  and  cold  gas  are  in  conflict.  A cold  gas  high  pressure  ejector 
in  this  speed  range  requires  a converging  (accelerating)  inlet  and  a diverging 
(decelerating)  outlet,  while  a hot  gas  low  pressure  ejector  requires  a diverging 
(decelerating)  inlet  and  a converging  (accelerating)  outlet.  Thus  the  ejector 
having  high  pressure  and  temperature  injected  gas  is  at  a conflicting  design 
point  where  the  best  performing  ejector  is  virtually  a constant  area  duct.  To 
counteract  this  difficulty  it  is  possible  to  utilize  ejector  design  criteria 
dictated  by  the  second  solution  with  subsonic  mixing  o£  to  use  a turbojet  with 
small  pressure  ratio  and  high  temperature  as  illustrated  on  Figure  9,  or  alter- 
natively to  use  a high  pressure  ratio  and  low  temperature  effluent  gas  of  the 
fzmjet  type. 
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d.  Rocket/Ejector 

The  use  of  an  ejector  in  combination  with  a rocket  type  gas  generator 
offers  performance  characteristics  similar  the  those  of  the  turbojet/ejector 
combination  as  shown  on  Figure  10.  Since  the  rocket  engine  does  not  require 
intake  of  ambient  air,  the  thrust  augmentation  (4>) , was  calculated  without 
ram  drag  of  the  primary  jet,  but  with  proper  consideration  of  the  ram  drag 
of  the  induced  flow.  Perfonn*mce  is  good  at  high  flight  speeds  under  all  design 
conditions,  particularly  when  the  ejector  is  designed  under  the  second  solution. 

As  indicated,  subsonic  mixing  provides  superior  performance  to  that  with  super- 
sonic mixing.  It  is  important  to  note  that  at  low  speeds,  the  thrust  augmentation 
of  this  system  is  extremely  large  provided  the  design  configuration  is  that  which 
corresponds  to  the  second  solution  with  subsonic  mixing. 
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SECTION  V 

REAL  FLUID  EJECTOR  EFFECTS 

The  realistic  evaluation  of  the  performance  of  an  ejector  thruster  depends 
upon  a knowledge  of  the  correct  values  of  the  parameters  which  describe  the 
deviation  of  the  flow  from  that  of  the  ideal,  completely  mixed  flow.  Certain 
of  these  parameters  can  be  evaluated  by  analytical  methods  using  conventional 
boundary  layer  theory  while  others  require  the  use  of  experimental  data.  The 
methods  available  for  evaluation  of  these  parameters  are  described  in  the 
following  discussion. 


1.  Nozzle  efficiency  (h^j) 

The  nozzle  efficiency  represents  a measure  of  the  loss  of  momentum  occuring 
in  the  primary  nozzle  duct  of  the  ejector.  This  loss  can  be  evaluated  using 
standard  duct  loss  factors  which  are  well  documented  for  certain  configurations 
such  as  corners,  diverging  and  converging  cross-sections  etc.  However  certain 
types  of  nozzles  which  are  utilized  in  ejector  design  involve  undocumented  con- 
figurations and  it  is  therefore  desirable  to  evaluate  the  loss  experimentally. 

To  accomplish  this  it  is  essential  to  operate  the  primary  jet  as  a free  stationary 
jet  and  to  measure  its  thrust,  mass  flow  and  stagnation  properties.  From  this 
experiment  and  the  ideal  thrust  (F.,  ,),  calculated  using  the  lossless  reservoir 

pressure 


ideal 


m 2c  T 

pV  p op 


[-fen 


where  F = the  stationary  thrust 

m - the  mass  flow  rate 
P 

T = the  stagnation  temperature 
op 

n = (7-1) A 

= reservoir  pressure 

the  nozzle  efficiency  can  be  evaluated  as 


N F 

ideal 
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When  supercritical  pressure  ratios  are  utilized,  the  ejector  requires  a nozzle 
which  is  converging  only  since  the  exp^msion  to  the  throat  pressure  can  occur  isen- 
tropically.  In  this  case  the  ideal  thrust  of  the  nozzle  as  a free  jet  is 
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^ideal 
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pV  P op 


b-iU] 


Measurement  of  the  thrust,  mass  flow  and  stagnation  temperature  then  can  be 

utilized  to  determine  ri  since 

N 


n = — 

^ F* 

ideal 


where  for  air. 


m /2c  ' 

pV  P 


= 0.5283 


If  the  primary  jet  is  a wall  jet,  additional  losses  can  occur  due  to  skin 
friction  which  can  be  evaluated  by  employing  boundary  layer  theory.  Since  the 
ejector  performance  <)>,  is  almost  directly  proportional  to  n^j,  using  wall  jets 
as  primary  jets  should  be  avoided. 
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2.  Inlet  loss 

The  Inlet  loss  of  the  induced  flow  occurs  at  the  entrance  to  the  ejector  as 
a result  of  blockage,  inlet  separation  and  skin  friction.  Since  in  any  real  ejector 
these  losses  cannot  in  general  be  evaluated  by  2uialytical  methods,  it  is  necessary 
to  utilize  an  experiment  in  which  the  distribution  of  pressures  and  temperature  are 
measured  at  a station  immediately  downstream  of  the  inlet.  Then  if  at  the  survey 
station 

p “ static  pressure 

=■  stagnation  pressure 
T^  = stagnation  temperature 
the  average  static  pressure  is 


p = ^ //  p dy  dz 


(79) 


the  entrainment  ratio  is 


dz 


and  the  queintity 
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With  the  )cnowledge  of  the  primary  jet  characteristics,  (m 

F 

n ))  and  the  ejector  geometry,  the  inlet  area  ratio  is. 
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The  value  of  can  then  be  determined  by  an  iterative  solution  of  the 
following  equations. 


(83) 


(84) 


(85) 


i • 


I 


-50- 


I 


2 C T ^ 

p o" 


b-m 


m 

3,  = - (v  + rU,  ) 
1 ' pi  1 ' 


+ p + PU 


(conservation  of  momentum) 


With  the  solution  for  p 
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is  also  determined  as  shovm  above  emd 
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as  described  in  Equation  8,  if  Ug/U^  = 1 

Experimental  determination  of  the  inlet  drag  coefficient  for  the  STAMP  AJDE 

(Small  Tactical  Aerial  Mobilitly  Platform,  Alperin  Jet-Diffuser  Ejector)  was 

performed  at  FDRC,  using  an  assumed  uniform  total  temperature  equal  to  the  mass 

flow  averaged  total  temperatxire  of  primary  and  induced  flows.  The  results  of 

this  experimental  analysis  yielded  inlet  drag  coefficients  (C^^) > ranging  from 

0,013  to  0.016,  for  the  jet-diffuser  ejector,  depending  upon  primary  nozzle  design 

and  attitude.  Since  the  AJDE  has  throat  dimensions  of  10.16  cm  x 38.1  cm,  the 

two-dimensional  value  of  C, . is  alDOUt  0.010  to  0.013. 

di 

At  supersonic  flight  speeds  the  inlet  loss  is  primarily  due-  to  inlet  shoc)c 
waves,  therefore  it  is  more  conveniently  expressed  as  the  inlet  recovery  factor 


’^i  = F 


This  factor  depends  upon  the  wave  configuration  of  the  inlet  cmd  its  max- 
imum acceptable  value  for  decelerating  inlets  of  the  type  used  for  engine  inlets 
is  specified  by  the  conventional  design  requirement  set  forth  in  MIL-E-5007D,  as 


'^i 


for  1 < M <5. 


= 1 - 0.075  (M„  - 1) 


t 


since  ejectors  translating  at  supersonic  speeds  may  have  subsonic  or  super- 
sonic mixing  (Mj^  may  be  greater  or  less  than  1.0),  the  inlet  may  be  subject  to  the 
shock  wave  loss  reflected  in  the  factor  (0^)  only  when  < 1.  If  the  ejector  is 
designed  to  utilize  supersonic  mixing  (Mj^  > 1) , the  inlet  may  be  either  deceler- 
ating or  accelerating  but  in  either  case  since  a normal  shock  can  be  avoided  the 
inlet  loss  is  primarily  comprised  of  skin  friction  2md  the  value  of  must  be 
determined  experimentally  as  described  above. 

The  influence  of  as  given  by  MIL-E-5007D  on  the  thrust  augmentation 
of  an  ejector  translating  at  a Mach  No.  of  2 with  a ramjet  type  gas  generator 
is  illustrated  on  Figure  11,  where  as  shown  the  performance  remains  high  despite 
the  large  loss  due  to  shock  waves  at  the  inlet.  Under  the  first  solution,  the 
use  of  subsonic  mixing  provides  superior  ideal  performance  to  that  of  supersonic 
mixing  but  shock  losses  required  for  subsonic  mixing  result  in  a reduction  of  the 
optimal  value  of  and  a degradation  of  performance  which  then  becomes  equivalent 
to  that  for  supersonic  mixing  (see  Figure  11) . 
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Thrust  Augmentation, 


3 . Flow  non-uniformity  due  to  incomplete  mixing 


Under  the  assumption  of  incompressibility  it  has  been  shown  (Reference  3) 
that  the  performance  of  an  ejector  depends  upon  a generalized  diffuser  area 
ratio  D (=  6 /(I  + U^  /U^) 1 which  includes  the  kinetic  non-uniformity  terms.  In 
that  case  the  thrust  augmentation  is  uniquely  related  to  the  pressure  at  the 
injection  station  (p^^)  and  the  flow  pattern  is  a function  of  D.  Thus  as  the  flow 

non-uniformity  increases  the  performeince  and  the  generalized  flow  pattern  can 

* 

remain  unchanged  if  the  effective  diffuser  area  ratio  (6  ) is  increased  by  the 
factor  (1  + • 

In  a compressible  flow  the  analysis  is  more  complex  as  a consequence  of  the 
existence  of  thermal  as  well  as  kinetic  mixing  in  addition  to  the  compressibility 
effects.  In  general  it  can  be  observed  however  that  in  a compressible  flow  the 
varying  degrees  of  mixing  or  flow  non-uniformity  can  influence  the  choking  limits 
and  the  performance  of  an  ejector.  These  effects  can  be  observed  as  follows. 
Writing  the  expressions  of  Equations  34  to  37  in  a more  convenient  form. 
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It  has  previously  been  shown  that  at  the  choking  limit  the  discrimincuit 
2 

(1  - 4AC/B  ) becomes  zero.  This  implies  that  the  non-uniformity  terms  in  Equation 

2 —2 

92  can  influence  the  extent  of  the  choked  region.  For  example,  since  /U^  and 
C are  always  positive  an  increase  in  their  magnitude  results  in  an  increase  in 

F 2 

4AC/B  or  a decrease  in  the  magnitude  of  the  discriminant.  Obviously  this  results 
in  a widening  of  the  choked  region  or  the  requirement  for  operating  the  ejector 
at  values  of  further  from  1.0.  Increasing  thermal  non-uniformity  being 
directly  related  to  the  magnitude  of  p2^2^P2^2'  ® decrease  in  the  magnitude 
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of  4AC/b‘’  or  an  increase  in  the  magnitude  of  the  discriminant  since  p^V^/p^U^  is 

generally  negative.  This  results  in  a narrowing  of  the  choked  region  and  the 

possibility  of  operating  the  ejector  at  values  of  closer  to  1.0.  Since  the 

term  B/2A  includes  the  same  non- uniformity  term  in  its  denominator  its  influence 

2 

on  ^2  oppKSsite  to  that  of  4AC/B  . 

To  examine  the  influence  of  flow  non-uniformity  on  ejector  thrust  augmen- 
tation, it  can  be  shown  that  from  Equation  22  and  31 
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Since  - 4AC/B^  is  generally  small  compared  to  1,  from  Eq.  C-5  of  Appendix  C 
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In  view  of  the  fact  that  (1  - 4AC/B^)  = 0(u:^/U?)  near  the  choking  limit, 

^ ^ 2—2 

the  last  two  terms  have  orders  of  magnitude  higher  than  U'  /U^  and  are  neglected. 


As  shown  in  Appendix  C,  using  the  above  approximation  and  Equation  92,  the 

1/Y 

expression  for  X2(P2/P„)  near  the  choking  limit  becomes 
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Examination  of  Equation  95  indicates  that  increases  of  non-uniformities 
decrease  the  term  (P2/Poo^  influence  on  X^d  + /U^  t ^3/p3^3^ 

is  diminished  by  the  non-uniformity  terms  at  Station  3 in  the  denominator  of 
Equation  96.  Thus  since  the  thrust  augmentation  is  directly  influenced  by  the  term 

2 —2  . I , 

X^d  + u'  /Uj  + Pj  U'/p^U^) , as  indicated  by  Equation  56,  the  net  effect  of  flow 
non-unformity  on  <t>  Is  small. 

From  Equation  52,  for  an  isentropic  ejector  outlet. 


(97) 


Examination  of  Equation  97  in  the  light  of  Equation  95  indicates  that 
increases  in  the  kinetic  non-uniformity  at  Station  2 results  in  a decreased  value 
of  the  effective  outlet  area  ratio  (6*)  . Since  the  right  hand  side  of  Equation  96 

is  virtually  independent  of  non-uniformity  as  discussed  above,  an  increase  of 

* 

kinetic  non-uniformity  at  Station  3 increases  the  required  value  of  6 due  to  a 
decrease  of  X^.  Thus  for  an  ejector  with  a decelerating  outlet,  the  non-uniformity 
term  at  Station  3 is  larger  than  that  at  Station  2.  This  requires  an  increased 

it 

value  of  6 with  increasing  kinetic  non-uniformity.  With  an  accelerating  outlet 
this  situation  is  reversed. 
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4.  Correlation  with  Experiment 

The  influence  of  skin  friction  and  incomplete  mixing  between  Stations  1 and 
2 are  represented  by  the  last  three  terms  in  the  coefficient  A of  the  solution  for 
as  indicated  by  Equation  35.  If  it  is  assumed  that  a coefficient  represents 
the  entire  influence  of  skin  friction  (which  contributes  to  entropy  production 
between  Stations  2 and  3)  and  incomplete  mixing  between  Stations  1 and  3,  then 
Equation  35  tcikes  the  form 


where , 


A = lii  + ic 
* 2y  2 S 


S = (V^2’ 


(98) 


= total  ejector  length  after  injection 

and,  as  a consequence  of  isentropic  flow  between  station  2 and  3 (s^  = s^) , 
Equation  48  can  be  expressed  as 

, 2 


and  from  Equation  (28) 


Poo  “ ^2  \ ^Dl  / \ ^ / Poo 


(99) 


(100) 


2 \ pi 

Using  this  approximation,  which  incorporates  the  influence  of  incomplete 
mixing  into  the  coefficient  C^,  and  applying  the  first  solution  to  an  ejector 
whose  geometry,  injected  flow  and  nozzle  characteristics  correspond  to  those 
quoted  in  Reference  4 for  the  ARL  Ejector  provides  excellent  agreement  between 
theory  and  Configuration  F performance  data  of  Reference  4,  with  = 0.006. 

This  high  performing  ejector  had  a diffuser  in  which  the  flow  remained  attached 

•k 

(i5  = 6)  . Configuration  D which  is  reported  to  have  had  the  Scune  length  after 
injection  as  that  of  Configuration  F but  which  had  a longer  mixing  section  euid 
a shorter,  wider  angle  diffuser  displayed  performance  which  was  considerably 
lower  than  that  of  Configuration  F and  also  lower  than  the  cdxsve  described 
analysis  (at  6 > 1.4)  as  shown  on  Figure  12.  Further  the  inception  of  diffuser 
separation  is  made  evident  by  the  drop  in  value  of  U^  with  increasing  geometric 

diffuser  area  ratio  shown  in  Reference  4.  It  can  therefore  be  concluded  that  the 

* 

Configuration  D diffuser  flow  separated  cind  that  6 <6. 
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Since  the  ARL  ejectors  had  a throat  dimension  of  10"  x 60",  can  be 
corrected  to  represent  the  coefficient  of  skin  friction  based  on  the  wetted 
surface  area.  This  value  is  0.0051  (=0.006  x 6/7),  which  is  very  close  to 

the  coefficient  of  skin  friction  of  a flat  plat  turbulent  boundary  layer  over 

£ 

a wide  range  of  Reynolds  numbers  near  10  . 

This  example  demonstrates  that  for  a solid  diffuser  ARL  ejector  with  length 

to  throat  ratio  mixing  of  the  primary  and  the  induced  flow  is 

adequate,  and  the  ejector  performance  is  controlled  by  factors  n , C , . , skin 

N di 

friction  (C^)  and  most  importantly,  the  effective  diffuser  area  ratio  6*. 


'i  < 


\ 
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SECTION  VI 


JET  DIFFUSER  EJECTORS 

Jet  diffusers  are  as  the  name  implies,  diffusers  which  are  comprised  of 
fluid  jet  sheets  which  can  support  a pressure  discontinuity  and  thereby  perform 
the  function  of  a solid  surface.  Thus  the  actual  length  of  the  solid  surfaces 
can  be  drastically  reduced  while  achieving  the  same  recovery  of  jet  kinetic 
energy  or  leaving  loss.  Recent  designs  of  jet  diffusers  by  FDRC  (Reference  1) , 
have  utilized  the  diffuser  jet  for  several  purposes  which  aid  in  the  improvement 
of  ejector  performance. 

By  injection  of  the  diffuser  jet  at  the  appropriate  location  within  the  solid 
diffuser  surface,  it  can  prevent  separation  despite  the  use  of  extremely  large 
divergence  angles.  Tests  have  demonstrated  that  with  a diffuser  jet  it  is  possible 
to  avoid  separation  with  solid  wall  divergence  as  large  as  60  degree  half-angle. 

In  addition,  when  the  diffuser  jet  emerges  from  the  end  of  the  solid  surface  with 
a positive  angle  (B)  (Figure  13) , and  with  a stagnation  pressure  or. 

than  it  produces  a jet  diffuser  which  increases  the  effective  diffuser  ai  a 
ratio  (a)  to  a value  which  is  larger  than  the  geometric  diffuser  area  ratio  -) . 

A further  indirect  advcintage  accruing  from  the  use  of  a diffuser  jet  is  related 
to  the  increased  effective  length  of  the  diffusion  process  which  thereby  permits 
the  mixing  of  primary  and  induced  flow  to  continue  effectively  downstream  of  the 
solid  surfaces. 

Jet  diffusion  has  been  shown  to  provide  the  above  described  advantages  in 
subsonic  flow  and  its  mass  flow  emd  momentum  requirements  for  optimal  performance 
have  been  determined  emalytically  (Reference  3)  and  experimentally  (Reference  1) . 
For  cui  ideal  jet  diffuser  ejector,  the  ejector  flow  will  be  fully  mixed  prior  to 
the  return  to  ambient  pressure  (far  down  stream  of  the  ejector  exit) , therefore, 
the  following  analysis  is  based  upon  the  assumption  that  the  flow  in  the  jet 
diffuser  ejector  is  uniform,  and  introduces  an  empirical  const£mt  to  correct  the 
imperfection  of  the  jet  diffuser. 

The  governing  equations  upstream  of  station  2 are  identical  to  those  for  a 
solid  diffuser  ejector,  except, 

1)  flow  is  assumed  uniform,  non-uniformity  terms  vanishes, 

2)  C_  = 0,  due  to  the  non-existence  of  a mixing  duct, 

3)  skin  friction  on  solid  surface  of  diffuser  is  accounted  for  in  the 
momentum  flux  of  the  diffuser  jet. 
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1.  Jet-Diffuser 


Previously,  in  the  discussion  of  the  flow  through  the  ejector  it  has  been 
assumed  that  the  flow  returned  to  ambient  pressure  at  Station  3,  which  location 
coincided  with  the  end  of  the  effective  diffusion  process  in  a solid  diffuser 
ejector.  If  however  a jet  diffuser  is  utilized,  the  flow  returns  to  ^unbient 
pressure  far  downstream  of  the  end  of  the  solid  surfaces.  Between  the  end  of  the 
solid  diffuser  and  the  region  where  the  flow  returns  to  ambient  pressure  there 
exists  a curved  jet  sheet  called  a diffuser  jet.  In  this  case  the  outlet  end  of 
the  ejector  appears  as  illustrated  on  Figure  13. 


@ (3)  ® 


Figure  13.  Jet-Diffuser 


The  diffuser  jet  is  injected  as  a wall  jet  at  Station  J,  and  encounters  a 
skin  friction  along  the  solid  surface  until  Station  3 where  the  solid  surface 
terminates  at  an  angle  6 with  respect  to  the  eucis  of  symmetry. 

The  momentum  theorem  between  Stations  3 and  4 may  be  written  as  follows. 


P3X3  + 2 / p_dy  + m^U,  + m^V^cos  B = p^X^  + m^  + ih^V 


X3/2 


c 3 d d 


■44  c 4 


d d 


(101) 


or  since  p^  » p^,  this  may  be  expressed  in  a dimensionless  form  as 


<P3  - PJ  2 ^ 2 


A V , 
P pi 


ro  V - 
P pi  X,/2 


*4/2  mV 

/ (p  - pjdy  + (1  + r)(X  -X  ) - (1  - cos  6)  = 0 (102) 

“3/2  ■ ^ p pi 
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The  Integral  involves  p_  (the  pressiire  on  the  external  streamline  of  the 
diffuser  jet) , which  equals  p^  in  the  stationary  case  and  in  that  situation  the 
integral  vemishes.  However,  in  the  translating  case  the  integral  ceui  be  evaluated 
by  replacing  the  jet  sheet  with  a vortex  sheet  as  is  done  in  jet  flap  theory. 

Thus, 


U„  - U* 


(103) 


U = U.  + U’ 

+ 4 

where  U*  is  assumed  small  compared  to  and  U^,  and  by  applying  Bernoulli's 
Equation  to  each  side  of  the  jet  sheet,  eind  assuming  identical  density  for  external 
and  internal  flows, 

P»  - p+  = PV' 

P_  - Poo  = PU’U„ 


and  therefore. 


P_  - = PU'  (U„  + U^) 


(104) 


where 


neglecting  second  order  terms  in  the  perturbation  velocity. 

The  momentum  law  in  the  y-direction  for  the  upper  isolated  jet  sheet  is 
00 

/(p_  - p^)dx  = p^tV^(sin  0-0) 

X 

t = the  jet  sheet  thickness 

0 = the  local  angle  w.r.t.  the  thrust  direction  (x) . 
Differentiation  with  respect  to  x yields 

.2 

= P_  - Pa 


(105) 


since 


and 


R ^ + 

R = Radius  of  curvature 

sin  0 = 


'1  + y 


d 3in0 
dx 


(1 


R 


(106) 
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Using  Equations  104  and  106 


P^tV^/R  = PU- (U„  + U^) 


and  since 


“d'^d  = 

= velocity  of  diffuser  jet 
= mass  flow  of  diffuser  jet  (2  sides) 


Therefore 


2P  RCU^  + UJ 


Using  Equation  110,  the  value  of  p_  - can  be  found  from  Equation  104  to  be 


P_  - Pco  = 


Vd^«> 

2R(U4  + UJ 


This  permits  integration  of  the  integral  in  Equation  102. 


/(p_  " pJ<iy  = 


d d °°  I ^ _ d d °° 
2(^4  + UJ  X3/2  ^ ^ U4  + 


(1  - cos  6) 


other  terms  in  Equation  102  are  expressed  conveniently  as  follows. 
Similar  to  Equation  (29)  the  p^  term  can  be  expressed  as 


^3^3  nJ  ,, 

1^77=  - ^ 

P pi 


^4 


and  noting  that 


m V , 
P pi 


Substituting  Equations  112,  113  and  114  into  Equation  102  results  in  the  ideal  jet 
diffuser  equation  which  relates  to  X^,  or. 


^ _ Y+l  \ . nJ  na6 

4 “ 2Y  3 2\^  ~ 2(l+r)  /P^  \n 


\Pl/  1 I" 

(P  \n  1+r  m V , 


(1  - cos  6) 


44:) 


. ...i. 


where  can  be  calculated  from  by  Equation  50,  if  6 is  known,  or 


3 \ * j y '2\ 

By  substituting  subsript  4 for  subscript  3 in  Equation  48  and  under 
the  uniform  flow  and  isentropic  outlet  assumption 

.2 


(116) 


(117) 


The  effective  diffuser  area  ratio  a can  be  obtained  by  assuming  isentropic  flow 
in  the  diffuser.  Substituting  0 for  6 emd  subscript  4 for  subscript  3 in 
Equation  50  neglecting  all  the  non-uniformity  terms,  and  noting  that  s^=  s^, 


results  in 


■ 4 


■ 4 


_1_ 

y-T 


(118) 


Combination  of  Equations  117  and  118  gives  the  relation  between  0 and  X^,  or. 


= ^2 


(r 


(119) 


Eliminating  X^  from  Equations  115  and  117  and  noting  that  all  other  terms  can 
be  expressed  in  terms  of  (p, /p..),  gives  the  solution  of  the  ideal  iet-dif fuser 
ejector  problem. 

In  the  realistic  situation,  imperfections  do  exist  in  the  design  of  jet- 
diffusers.  For  example,  in  a three  dimensional  ejector,  the  peripheral  variation 
of  the  diffuser  jet  thickness,  and  the  non-uniform  peripheral  distribution  of  3, 
Ccui  degrade  the  performance  of  a jet-diffuser  and  result  in  a decrease  of  the 
effective  diffuser  area  ratio  (0)  and  a premature  termination  of  the  effective 
mixing  process  which  then  decreases  the  entrainment  capability  of  the  ejector, 
euid  penalizes  the  ejector  perfomumce.  Also  the  definition  of  outlet  area  ratio 
5 (=  X^/Xj)  is  arbitrary  since  the  flow  is  highly  two-dimensional  (90®  included 
angle  in  the  AJOE)  it  is  appropriate  to  introduce  an  empirical  constant  to  relate 
the  ejector  geometry  to  the  outlet  area  ratio. 
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The  ii««)erf actions,  attributable  to  the  peripheral  non-uniformity 
dimensionality  and  three-dimensionality  of  the  flow  due  to  finite  lon( 
dimensions  can  be  described  by  an  efficiency  factor  n in  the  form. 


or  replacing  6 with  6 n 


(122) 


The  complete  solution  of  a jet-diffuser  ejector  can  be  obt 
from  equations  117  ^md  121  and  solving  for  (p^^/p^)  if  6 
for  6 if  (p./p^)  is  given. 


(123) 


xne  ratio  or  airruser  jet  to  primary  jet  momentum  which  occurs  in  Equations  123 
and  121  depends  upon  the  degradation  which  occurs  as  a result  of  the  skin  friction 
encountered  by  the  diffuser  jet  along  the  solid  surface  between  station  J and 
Station  3.  This  skin  friction  force  can  be  evaluated  using  the  Blasius  Law 
when  the  Reynolds  number  of  the  flow  is  less  than  5 x 10^  (laminar  flow)  or  by 
using  the^Prandtl-Schlichting  Law  for  Reynolds  numbers  greater  than  or  equal 
to  5 X 10  (transition  to  turbulent  flow) . 


The  Reynolds  number  based  upon  length  of  wetted  surface  is  (Reference 


Re/Jl,,= 


^'^da>  ^od”d» 


VW-DV 


od 


^ <T^/T  ^)  + (198.6/T  ,) 

y-l  <3°°  od od 

1 + (198. 6A  j) 
od 


M^=  12/ (Y  - 1)H(P^/Pj"  - 1] 


3/2  -8  Ibf  - sec 

y ^ = 2.270[T  , ^ /(T  , + 198.6)1  ^ 10 
od  od  od 


2 2 o 

c = 4290  ft  /sec  R 

V 


ft 


To  determine  we  used  the  geometry  of  a family  of  jet  diffuser 
aD 

ejector  designed  by  FDRC  whose  geometry  is  described  by. 


^dj'^^2 


6 < 1.0304 


= (1.09996  - 1.1333)  1.0304  < S < 1.5855 

d3  2 

Z^./X^  = (0.70716  - 0.5105)  6 > 1.5855 

d3  2 — 

= 4 in.  = 10.16  cm. 

Thus  if  is  in  psia,  X^,  in  inches,  and  T^  in  °R  then  for  air 

Re^  = 1.51  X 10^M^P^/T^)(!l^.A2)(T^/T^)^*^[(T^yT^^)  t (198.6/T^^)  ]X. 

Then,  according  to  the  Blasius  Law  (Reference  6) 

Cfdj  - 1.328/Y^  Re^  < 5 X 10^ 

and  according  to  Prandtl-Schlichting,  (Reference  6) 

7 CO  5 

Cfdj  = 0.455/(log  Rej^)  - (1613, S/Re^^)  Rejj^  > 5 x 10 


and  the  compressibility  correction  term  is  (Reference  6) 

{1  + 0.85(M^)((Y  - l)/2]} 
fd3 , comp  fd j , incomp  d» 

Using  Cj^j  determined  in  this  manner  permits  evaluation  of 

m.VVm  V , =>  (m.V,  /A  V ,)  [1  - <<^.,u/sJC 
d d'  p pi  d d«f  p pi  dj  “ fdj 


= (p^Pj^)  (sya^) 


II  - I^dj/"oo)Cfdjl 


(124) 

(125) 

(126) 

(127) 

(128) 


(129) 

(130) 

(131) 
since 

(132) 
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2.  Correlation  with  Experiment 

To  illustrate  the  validity  of  this  analysis,  it  was  used  to  evaluate  the 
performemce  of  a series  of  jet-diffuser  ejectors  which  were  tested  at  FDRC 
under  the  STAMP  (Small  Tactical  Aerial  Mobility  Platform)  Program  (Reference  1) . 
These  ejectors  were  three-dimensional,  having  widely  diverging  diffuser  surfaces 
at  the  sides  and  flat,  non-diverging  end  surfaces.  The  lack  of  divergence  of  the 
ends  of  the  diffuser  resulted  in  a collapse  of  the  flow  since  the  non-diverging 
diffuser  jets  could  not  sustain  the  pressure  gradient.  To  delay  this  collapse, 
end  plates  of  various  sizes  were  added  to  extend  the  solid  surfaces  beyond  the 
solid  surfaces  at  the  sides  of  the  ejector  as  illustrated  on  Figure  14.  Measure- 
ments of  thrust  augmentation  were  made  for  this  series  of  end  plate  configurations 
with  other  geometric  and  thermodynamic  parameters  held  constant.  This  data  is  also 
presented  on  Figure  14. 

Using  the  geometric  and  thermodynamic  parameters  corresponding  to  those  of 
the  STAMP  ejector  tests,  the  values  of  and  were  calculated  using  the 

above  method.  The  improvement  of  jet  diffuser  efficiency  (^^j)  with  increasing 
size  of  the  end  plate  is  illustrated  on  Figure  15. 
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SUCTION  VII 


CONCLUSIONS  AND  R1':MARKS 


Tho  poi'fonnaiico  of  systoms  wliich  propel  airborne  vehicles  can  bo  qroatly 
improved  by  tho  use  of  combination  thrusters  comprised  of  qas  generators  and 
eiectors,  at  any  present  doy  practical  operational  condition.  Tho  maximum 
advantage  derivable  from  these  combination  thrusters  can  only  bo  achieved  as 
the  result  of  a careful  optimization  of  tJio  characteristics  of  the  gas  generator 
and  the  ejector  in  view  of  tho  operational  requirements. 

At  low  speeds,  optimal  overall  perfontivinco  is  achieved  with  gas  generators 
whicl>  supply  relatively  cold,  high  pressure  ratio  gas  to  tlie  ejector.  Such  gases 
can  be  obtained  for  ex/implo  from  tho  fan  air  of  a turbofan  engine.  High  e joct.or 
performance  can  also  be  achieved  with  hot  gas  injection  at  low  pressvire  ratios 
but  the  overall  performance  of  .such  systems  may  suffer  from  tlie  poor  thennal 
efficiency  of  the  gas  generator  at  low  flight  speeds. 

At  high  subsot\ic  or  sujiorsonic  speeds,  high  overall  performance  can  be 
achieved  by  comliination  thrusters  which  suijply  very  hot  injected  gas  to  tlie 
ejector.  This  is  a vefi'  practical  arrangement  since  at  those  high  speeds,  a 
' gas  generator  can  supply  hot  effluent  gas  with  high  thermal  efficiency  and  with 

an  extremely  simple  and  reliable  gas  generator  of  the  ramjet  type. 

The  performance  degradation  resulting  from  nozzle  duct  loss,  inlet  drag, 
skin  friction  and  incomplete  thermal  and  kinetic  mixing  has  been  considered  in 
the  analysis  and  methods  have  boon  discussed  for  evaluation  of  the  loss  coef- 
ficients. Tho  performance  of  an  ejector,  after  consideration  of  the  losses,  has 
tH<,>n  shown  to  agree  very  closely  with  experiments  performed  at  ARl.  and  I'DKC 
uiuler  stationary  conditions. 

Since  no  data  are  available  for  translating  ejectors,  the  relative  importance 
of  each  loss  coefficient  during  translational  motion  of  tho  ojoctor  cannot  be 
ev.iliMteil.  liowevor,  a complete  treatment  of  tho  constant  area  mixing  ejector  was 
pie.K  ted  witli  all  loss  coefficients  for  future  design  and  analysis  of  ejector 

] 1 T '!■  l ems  . 

j * toi.»l  of  nine  types  of  ejector  configurations  are  required  to  accomodate 

1.  ':■>*»  patterns  wlilch  exist  for  optimal  ojoctor  performance  tuider  all  conditions 

been  examined.  Only  one  configuration  has  Inien  developed  to  date  and 
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therefore  only  this  one  configuration  has  been  used  for  comparison  of  theory  and 
experiment.  This  configuration  involves  an  acclerating,  converging,  subsonic 
inlet  and  a decelerating,  diverging,  subsonic  outlet.  Evaluation  of  the  loss 
coefficients  for  other  types  of  ejector  configurations  would  require  design, 
fedirication  and  testing  of  each  type  of  ejector. 

Using  the  design  principle  of  the  first  solution,  it  has  been  shown  that 
ejectors  cannot  operate  effectively,  in  the  medium  subsonic  speed  range,  when 
the  injected  gas  contains  both  thermal  and  kinetic  energy  since  the  optimal 
ejector  configurations  for  hot  gas  injection  are  very  different  from  those  for 
cold  gas.  Therefore  it  is  essential  that,  in  case  of  a turbofan  engine,  the 
efflux  from  the  fan  and  the  hot  gas  be  injected  into  separate  ejectors  to 
achieve  the  best  result. 

As  a result  of  the  analysis  based  on  the  conservation  laws  of  mass,  momentum 
and  energy  within  the  Second  Law  of  Thermodynamics,  the  use  of  ejectors  in  combi- 
nation with  gas  generator  engines  has  been  shown  to  provide  a means  for  producing 
more  mechanical  energy  than  that  which  is  injected  from  its  gas  generator.  This 
can  only  result  from  the  conversion  of  injected  thermal  energy  to  mechanical 
energy  in  the  ejector,  in  an  amount  in  excess  of  the  intrinsic  loss  resulting  from 
the  mixing  process.  Thus  the  otherwise  wasted  thermal  energy  of  a free  jet  can 
be  converted  to  mechanical  energy  similar  to  the  Brayton  cycle  process.  The 
conversion  of  lost  mechanical  energy  (leaving  loss)  to  useful  energy  has  also 
been  described  in  terms  of  the  relative  values  of  the  propulsive  efficiency  of 
the  free  jet  from  the  gas  generator  and  the  propulsive  efficiency  of  the  combi- 
nation thruster. 

It  has  been  shown  in  References  2 and  3 that  boundary  layer  ingestion  has 
an  effect  on  performance  similar  to  that  of  a reduction  of  ram  drag,  and  since 
<(>  is  a monotonic  function  of  q/AP,  in  an  incompressible  fluid,  a decrease  of 
ingestion  velocity  results  in  an  increase  of  thrust  augmentation.  Therefore  in 
an  incompressible  fluid  it  is  always  desirable  to  locate  the  ejector  at  a region 
where  boundary  layer  or  wake  can  be  ingested.  In  a compressible  fluid,  it  has 
been  shown  that  the  performance  dept  upon  temperature  and  pressure  of  the 


energized  fluid  and  does  not  always  decrease  with  increasing  flight  Mach  No. . 
Therefore  ingestion  of  boundary  layer  or  wake  fluid  while  decreasing  the  ram 
drag  also  decreases  the  inlet  compression  which  is  required  for  high  ejector 
perfona«mce  under  certain  conditions.  Therefore  optimal  ejector  integration 
may  require  full  usage  of  boundary  layer  ingestion  or  free  strejun  ingestion. 
In  the  latter  case  the  ejector  may  be  better  integrated  upstreeun  of  the 
boundary  layer  or  wcdce. 
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APPENDIX  A 


NORMAL  SHOCK  RELATION  BETWEEN  THE  FIRST  AND  THE  SECOND  SOLUTION 


In  the  ideal  case  when 
Equations  34  to  37 


0 and  flow  is  uniform  at  Station  2,  then,  from 


where,  A 


From  Equation  A-1 


Substituting  Equation  A-10  into  A-9 


Using  Equation  A-2  and  A-4  to  express  C and  A in  terms  of  J and  y 


Multiply  the  numerator  and  the  denominator  by 


If  the  first  solution  has  a M_  of  M_  and  the  second  solution  has  a M, 


which  is  a normal  shock  relation  (Reference  5) 


APPENDIX  B 


ENTROPY  EXPRESSION  FOR  A NON-UNIFORM  FLOW  FIELD 


Assumptions : 

1.  uniform  pressure 

2.  density;  p = p + p 


The  total  entropy  relative  to  Its  freestre2un  value  is 


then  Equation  B-3  becomes 


Neglecting  the  terms  with  order  of  magnitude  higher  than  2,  and  integrate 
Equation  B-5  becomes. 


= f fii  4. 


APPENDIX  C 


LINEARIZATION  OF  X 


NEAR  THE  CHOKING  LIMIT 


Then 


(r-l)  (2T-1)  Or-l)  ^4,  4^  ...1 

o 24  J 

(1  i X)  { [i  + X - T x'-  (T.zi).L2y-i)  (.3y-i)  ,4^  . . .] 

+ y[Tx-  (y-l)x^  ? ^3_  (lr.lH2X-lJ13Y-.l)  ^4^  ^ ...| 

[|l  _ lli  T iliimziL  (Y^-1)  (Y-1)  (2Y-1)  ^^4,  ...J 


+ y[+  X - Yx^  •••]+  •••  

Since  near  the  choking  limit,  O(x^)  = 0{y)  = O^U^^/U^j, 
therefore  the  leading  term  of  the  above  equation  is, 

(1  + x)[l  + Yx(l  + = 1 - ^ x^ 


= ill  + 2 + 2 


\ P2  / 

\ fi  + 5 ( izx  !!il  + SS  + 1 c \1 

7 (1-Y)/2  + z(Y+l)/2y  2 -2  ^ ^ 2 "'F^J 


(C-5) 
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2 / 2 — 2\ 

Near  the  choking  limit,  x = ^1^2  ^2)'  from  Equation  C-4, 
o(o'2/u2) 


Z = 1 + 


{C-7) 


therefore, 


= 1 + 


(1-Y)/2  + z(Y+l)/2 
and  Equation  C-6  becomes, 

+ X)  [1  + Yx(l  + y)] 


(C-8) 


(1 


1/Y 
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(C-9) 


Therefore, 


,2 

2 


p'  u' 

+ _£ =.  - — c 

2 ^ - 2 Si 


P U 
2 2 


J(1  + r)^(Y^  - 1)  (Y 


,2 


(l-Y  J(1  + r)^(Y^  - 1)  (Y  + 1)  1 /,  3-Y  ^2  , ^2  ^2  Y \ 

L op  pi  J 


(C-10) 


2 (T  /T  , 
op  pi 


and  Equation  (94)  becomes  (near  the  choking  limit) 

\t1/Yi 


xj^r  = 


no 


Ap, 


2(T  /T  - 1) 

SB—El 


-] 


(Y+1) /Y 


+ r/  (Y  + 1) 


I 


1-Y  ^ J(1  t r)^(Y^  - 1) (Y 


jLl)  Wi  . Izi  ^ + Ili  Cj\ 


2Y 


1 + rX,  + 


nOi 


1 2(T  /T  , 

op'  pi 


(C-11) 
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